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SUMMARY 


A laboratory experiment has been conducted in which acoustic signals propagating 
over both anechoic and reflective surfaces were detected by a pair of three-element 
microphone arrays. The phase component of the cross power spectrum between array 
elements has been examined for both ground impedance cases with the source location 
inferred from the slope of the phase spectra. The effect of ground reflections on 
both the accuracy and precision of position estimates made in this way is examined. 


INTRODUCTION 

The use of sensor arrays for the detection and localization of an acoustic- 
emitting source has found application in several fields. Sensors mounted on the 
surface of a material under stress have been used to detect acoustic emissions from a 
defect within the material (ref. 1). A multimicrophone acoustic ranging technique 
has been used to assess the degree of concentration of noise generated by the 
rotating blades of a wind turbine (ref. 2). Cross-correlation techniques between 
several pairs of microphones were employed to measure the sound radiated by a 
wing/flap configuration in a wind tunnel (ref. 3). Various applications generally 
make use of either correlation techniques for estimating time delay between array 
elements or phase information for the purpose of determining the direction of a 
coherent wave front as it passes the array. In general, the goal of either method 
has been to minimize the random error in the intersensor time delay (or phase differ- 
ence) in the presence of statistically independent noise at the array locations. It 
is usually assumed that dispersion of the estimates About their mean is due only to 
statistical sampling errors. In most applications, however, errors occur because of 
factors other than sampling considerations, such as wind and thermal refraction, 
scattering by objects or atmospheric inhomogeneities, and reflections from 
boundaries. 

The present report describes a simple approach used to locate an acoustic source 
passively and presents some accompanying experimental data to demonstrate its use 
under optimal conditions. The approach makes use of a linear best-fit regression to 
the slope of the phase-frequency plots between array elements. It was shown in 
reference 4 that phase data regression lines yield results with the same accuracy as 
other optimal time-delay estimation procedures. Moreover, phase data have an advan- 
tage in that extraneous noise sources having a different frequency content can be 
distinguished, and systematic errors can be assessed from the knowledge that the 
phase should be 0 for infinite wavelengths. 

The present analysis is for the case of two three-element microphone arrays. 
Tests to demonstrate the technique were performed in an anechoic facility in order to 
control meteorological and, hence, propagation effects. The objective was to estab- 
lish statistical bounds on range, elevation, and azimuth for a stationary source 
under optimum test conditions. Results are presented for both pure-tone and broad- 
band noise. Several source positions were tested in order to assess the effect of 
elevation angle on source localization. The effect of simulating extreme changes in 
boundary impedance on the ability to localize a source in space confidently is also 
addressed. 



SYMBOLS 


Values are usually given in SI Units but, where considered useful or expedient, 
are also given in U.S. Customary Units. Measurements and calculations were made in 
U.S. Customary Units. 


speed of sound, m/sec 

speed of sound at temperature T , m/sec 

phase speed in direction of line joining microphone positions x and 
y, m/sec 

acoustic frequency, Hz 
wave number vector, m " -1 

component of wave number vector in direction defined by microphone 
positions x and y, m 1 

slope of phase-frequency function (where phase is measured between 
microphones at positions x and y) 

range (distance from source to reference microphone), m 

distance from source to reference microphone in array I, m (ft) 

distance from source to reference microphone in array II, m (ft) 

separation distance between microphone arrays, m 

ambient temperature, K 

reference temperature, K 

source elevation angle, relative to plane of microphone array, deg 
elevation angle of source from array I, deg 
elevation angle of source from array II, deg 
microphone separation, m 

phase difference between signals measured at microphone positions x and y 
source azimuthal angle with respect to microphone array, deg 
azimuthal angle to source from array I, deg 
azimuthal angle to source from array II, deg 
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THEORY 


Assume that an acoustic wave of frequency f is propagating from a source 
located at spherical coordinates r , 0, and relative to the origin of a 

coordinate system, where r is the distance to the source, 0 is an elevation 
angle, and is an azimuthal angle, (See fig, 1,) 


Assume further that a three-element microphone array is located in the 0 = 0° 
plane of this frame of reference, with one microphone (microphone B) located at the 
origin and the other two (microphones A and C) located at r = £, xj) = 0° and at 
r = £, xp = tt/ 2, as shown in figure 2. 


If £/r << 1, the wave front curvature may be ignored. The component of the 
wave number vector in the 0=0° plane can be resolved into two orthogonal compo- 
nents, k and k , which correspond to the = 0° and xp = tt/2 directions, 
respectively. The difference in phase between signals received at microphones B and 
> BA ' is D ust for < 2 it. Similarly, c|> D „ = k_£ for k„£ < 2 tt. The 

>nas< 


A ' 


k BA 5 f ° r k BA 5 < 2ir * similarl y' ♦bc = k BC 5 for < 2tt. 

two pKase differences can be expressed in terms of frequency as follows: 


* -21SL 

BA °BA 


f =2^ 

50 C BC 


(la) 


(1b) 


where c BA and c B ^ are phase speeds in the 0=0°, = tt/2 and 0=0°, ^ = 0° 

directions , respectively. 

Let m BA and m BC be the derivatives of <j> and <f> BC , respectively, with 
respect to frequency. Then, 


m- 


BA 


= _2tt§ = ^BA 
"BA f 
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BC 
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Neglecting effects such as refraction, which may cause the ray paths to bend, the 
azimuthal angle ij; is 


= tan 



(3) 
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or 


ip = tan 


m 


BC 


in 


BA 


(4) 


and the elevation angle 0 is 


0 = cos 


,2 , 2 
-1 ( k BA + V 


From equations (2), 


m__f 

BA 


BA “ £ 


nif 

BC 


BC “ £ 


(5) 


(6a) 


(6b) 


Inserting equations (6) into (5) gives 


0 = cos 


-'U 


2 2 

m n , + ni 
BA BC 


(7) 


or 


0 = cos 


-1 c / 
2ttS|T 


2 2 

in + hi 
BA BC 


( 8 ) 


where c, the speed of sound, depends on ambient temperature T as follows: 


c = 



(9) 


where c Q is the speed of sound at temperature T Q . Insert equation (9) into (8) to 
express the temperature dependence of 0: 
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( 10 ) 



The bearing angles 0 and must be known at two different locations to 

determine the range to the source by triangulation. Assume a geometry as in figure 3 
in which two triangular arrays similar to that illustrated in figure 2 are located in 
the 0=0° plane. Array I contains microphones A, B, and C and array II contains 
microphones D, E, and F. Let the subscripts I and II refer to measurements made 
with respect to arrays I and II, respectively, and let s be the separation distance 
between microphones B and E, as shown in figure 3. Then, 


s sin ij; 

— — ■ — ■ ■ - '■ ' 

I cos 0^ sin( ^ 


(11a) 


and 


s sin i|^ 

II cos 0 sin( ^ + ifc^) 


(11b) 


Equations (4), (10), and (11) can be used to predict the location of a noise 
source from an analysis of relative phase spectra. Equations (1) indicate that the 
phase difference between signals received at two adjacent microphones in the array of 
figure 2 is directly proportional to frequency. Regression analysis on the phase 
spectra for microphone A relative to B and for microphone C relative to B will yield 
the slope values m BA and m B £, respectively, used to calculate bearing angles via 
equations (4) and (10). Bearing angles computed in the same way for a second array 
can be used to calculate range values as in equations (11). 


TEST DESCRIPTION 
Test Setup 

The microphone array geometry of figure 3 was established in the Langley 
Anechoic Noise Facility (fig. 4). Each array consisted of three 1/2-in. condenser 
microphones defining an isosceles right triangle, with the microphones that defined 
the equal legs of the triangle separated by a distance of 0.1 m (fig. 5). The micro- 
phone at the 90° corner of each array was defined as the reference microphone in that 
array, and all source locations were measured with respect to the center of the 
reference microphone diaphragm. 

The plane of each array was located 1 .79 m above the floor of the anechoic cham- 
ber. Covering the floor were polyurethane foam wedges 0.91 m in height, so that the 
plane of each microphone array was 0.88 m above the tips of the wedges. The spacing 
between the tips of the wedges varied from 0.3 to 0.4 m. An anechoic wedge treatment 
also covered the walls and ceiling of the chamber. The two arrays were separated by 
a distance of 5.26 m between reference microphones. Each array was located approxi- 
mately 1 .5 m from the sidewall wedge tips. 
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A small acoustic driver, fitted with a tapered flange, was used as a point 
source of noise (fig. 6). Four source positions were selected, as shown in 
table I. These source positions were determined relative to the reference microphone 
in each array by using a theodolite. The theodolite measurements are assumed to 
represent the true source location. 


TABLE I. - NOISE SOURCE POSITIONS 



The microphones were removed from their mounting hardware during theodolite measure- 
ments, and a specially designed theodolite mount was used to ensure that all measure- 
ments were made with respect to the position of the reference microphone diaphragm 
(fig. 7). 


Test Conditions 

Ten discrete frequencies were broadcast at each of the four noise source posi- 
tions in table I. The frequencies ranged from 200 to 2000 Hz in 200-Hz increments, 
with the lower frequency dictated by the performance of the chamber (not anechoic at 
frequencies substantially below 200 Hz ) and the upper frequency based on the micro- 
phone array geometry. (Wavelengths were chosen to be larger than the microphone 
spacing to ensure coherence.) In addition to the 10 pure tones, white noise in the 
range from 200 to 2000 Hz was broadcast at each source position. For each combina- 
tion of source position and frequency, data were acquired for a "soft floor" and a 
"hard floor" to simulate data acquisition under significantly different ground imped- 
ance conditions. The soft floor was achieved by leaving the anechoic foam rubber 
floor wedges in place, as shown in figure 4. For the hard- floor measurements, the 
floor wedges were removed, exposing a rigid plywood surface. 
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Data Acquisition 


Prior to the primary data -acquisition activity, a preliminary set of measure- 
ments was made to determine the phase characteristics of each microphone. The phase 
response of each microphone was compared with a reference microphone at each of the 
discrete-tone test frequencies. Differences in these phase responses for each micro- 
phone pair, which were averaged over the 10 pure-tone test frequencies, are presented 
in table II. These phase differences are small and were ignored in the test. 


TABLE II.- AVERAGE PHASE DIFFERENCE 
FOR EACH MICROPHONE PAIR 


Microphone 

pair 

Average phase 
deviation, 
deg 

Standard 

deviation, 

deg 

B,A 

0.3 

±0.5 

B,C 

.4 

+ .4 

E,D 

.7 

+ .8 

E,F 

1.1 

±7.0 

i . ... 


During the pure-tone measurements, a narrow-band analyzer was used to measure 
the amplitude and phase components of the cross power spectra for all three unique 
combinations of microphone pairs in each array. A sample rate of 8192 Hz and a band- 
width of 4000 Hz were used in these measurements. The amplitude spectra were dis- 
carded because only phase information is required to determine source position by 
the method of this paper. For the pure-tone measurements, only the one point in the 
phase spectrum corresponding to measurement frequency was meaningful, and the phase 
at that frequency was recorded for subsequent analysis. The broadband data were 
acquired in a manner similar to that of the pure-tone data, except that the entire 
phase spectrum in the range from 200 to 2000 Hz was recorded for subsequent analysis. 


RESULTS AND DISCUSSION 

For each combination of source position and ground impedance (four positions by 
two impedance conditions, resulting in eight cases altogether), phase differences 
were plotted against frequency for microphone pairs B, A, B,C, E,D, and E, F. This 
process was repeated for both the pure-tone and the broadband cases. Figures 8 and 9 
are the phase spectra corresponding to the soft-floor measurements of source posi- 
tion 1 (see table I) for the pure-tone and broadband measurements, respectively, 
whereas figures 10 and 11 are the corresponding phase spectra for the hard-floor 
case. These spectra are representative of the phase data acquired at the other three 
source positions. The hard-floor spectra deviate more from the theoretical straight 
line predicted by equations (1) than the soft-floor data, a result that is believed 
to be due to the influence of the hard-floor reflections on the phase measurements. 
Such reflections would cause signal cancellations and additions, depending on fre- 
quency, with a corresponding distortion in the phase spectra as seen in the data. 
Despite some fluctuations, all eight combinations of source position and floor imped- 
ance gave phase spectra that were well represented by a linear fit, as represented by 
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the straight lines in figures 10 and 11. Linear-regression correlation coefficients 
for the pure- tone data/ averaged over the four microphone pairs, are displayed in 
figure 12 for each combination of source position and ground impedance. The average 
correlation coefficients ranged from 0.83 to 0.96 for the hard-floor data and from 
0.92 to 0.99 for the soft-floor data, where a value of 1.00 would correspond to the 
theoretical prediction of equations (1). These results further suggest that source 
location errors are greater in the presence of reflections from a hard surface than 
from a soft surface. 


Accuracy of Acoustic Measurements 

The slope of the best straight-line fit for each phase spectrum was used in 
equations (4) and (10) to compute bearing angles to the source from the corner micro- 
phone of each array. These bearing angles were then used in equations (11) to 
compute the distance to the source. Figures 13 and 14 compare the hard-floor and 
soft-floor measurements of range, elevation angle, and azimuthal angle for arrays I 
and II, obtained from the pure-tone acoustic data, with the corresponding coordinates 
of table I, obtained by theodolite measurements. Figures 15 and 16 make similar 
comparisons between the broadband acoustic measurements and the theodolite 
measurements. 

The following observations can be made from the data in figures 13 to 16. The 
agreement between acoustic and theodolite measurements of range and bearing angles is 
better for the soft-floor case than for the hard-floor case. It is suspected that 
ground reflections contaminated the phase measurements that were used to forecast the 
source location. There was better agreement between acoustic data and theodolite 
measurements for the broadband data than for the pure-tone data, a result attributed 
to the more complete statistics resulting from the greater number of data points used 
in the broadband regression analysis (205 for the broadband data against 10 for the 
pure-tone data). The extent to which the broadband data improved the estimate of 
source position was greater for the hard-floor data than for the soft-floor data, and 
it was greater for the two lower elevation angles than for the two higher elevation 
angles. This again indicates the effects on phase measurements that hard-floor 
ground reflections had when compared with the reduced soft-floor ground reflections. 

The results of figures 13 to 16 are condensed in figure 17, which compares 
source locations measured from the acoustic data with source locations measured by 
theodolite. The comparisons are made in terms of the distance between source loca- 
tions measured acoustically and by theodolite. Figure 17(a) presents the results for 
the pure-tone data, and figure 17(b) presents the results for the broadband data. 
Again, there is better agreement with soft-floor measurements than with hard-floor 
measurements. The theodolite measurements agree better with the broadband acoustic 
results than with the sparser pure-tone data, with the extent of improvement greater 
for the hard-floor cases than for the soft-floor cases. 


Precision of Acoustic Measurements 

The upper and lower limits on the 90-percent confidence intervals for each slope 
estimate were computed during the regression analysis that fit the best straight line 
to the phase spectra of figures 8 and 9. Similar confidence intervals were computed 
for the other source positions. Equations (4), (10), and (11) were then used to 

transform these slope confidence intervals into confidence intervals in r, 0, 
and tp. The upper and lower limits on these spherical -coordinate confidence inter- 
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vals were converted to rectangular coordinates and plotted in figure 18 relative to 
the source position determined from theodolite measurements. Figures 18(a) and (b) 
correspond to one of the source locations and, for each figure, the 90-percent confi- 
dence intervals are represented as rectangles that define the 90-percent confidence 
regions in rectangular coordinate space. These 90-percent confidence regions are 
plotted for both the hard-floor and soft-floor data. Each figure illustrates both 
the horizontal plane and the vertical plane of the 90-percent confidence region. 

Also plotted in these figures are the hard-floor and soft-floor estimates of the 
source position. 

The magnitude of source position errors has been discussed earlier in this 
report. Figure 18 illustrates the direction of these errors as well as their magni- 
tudes. Figure 18 also graphically represents the precision with which source posi- 
tion estimates were made in these tests; it can be said with 90-percent confidence 
that the actual source position lies within the regions depicted in these figures. 

Note that the 90-percent confidence regions are smaller for the soft-floor data 
than for the hard-floor data, again indicating the effect of ground reflections on 
the precision with which an acoustic source can be located by the phase-slope 
method. The degree of reduction in the size of the confidence intervals going from 
the hard floor to the soft floor is greater at lower elevation angles than at higher 
elevation angles, a fact that is attributed to the greater effect of ground impedance 
and lower elevation angles. 

The 90-percent confidence intervals tend to be wider on the elevation angles 
than on the azimuthal angles, as indicated by the somewhat vertically elongated 
shapes of the confidence regions in figure 18. This result is attributed to the 
rather strong vertical temperature gradient that existed in the anechoic chamber 
during the noise measurements (typically 0.3°C per meter). This gradient is likely 
to have caused refraction effects that may have adversely affected the elevation 
angle measurements. A strong vertical temperature gradient would tend to bias the 
elevation angle measurements higher, a phenomenon that would be offset to some extent 
by the effect of ground reflections, which would bias the elevation angle estimates 
lower. Evidence of the temperature effect can be seen in figure 18 in which the 
90-percent confidence intervals in the vertical direction are generally larger for 
the higher elevation angles (for which the temperature variations are greatest) than 
for the lower elevation angles (for which the temperature variations are less over 
the propagation path). 

Bearing angles can be computed on the basis of phase information obtained at a 
single frequency, but the phase-slope method described in this paper offers signifi- 
cant advantages over such single-point measurements because phase information 
obtained simultaneously at several frequencies is used in the computations. Random 
distortions in the phase spectrum, caused by such factors as meteorological effects, 
reflections, and ambient noise, are manifested as positive and negative excursions 
about the best straight-line fit to the phase spectrum. Even spectra with relatively 
large fluctuations about the theoretical straight-line fit predicted by equations (1) 
can forecast bearing angles well, as long as such fluctuations are random in nature 
and do not systematically bias the slope of the phase spectrum in one direction or 
another, as illustrated by figure 11. These phase spectra were acquired under the 
relatively unfavorable condition of a high-impedance (rigid plywood) reflecting 
floor, with ground reflections interfering with the primary source signals to produce 
the distortions seen in these phase spectra. Since these distortions are more or 
less randomly distributed about the best straight-line fit, the corresponding regres- 
sion values still result in relatively good bearing angle forecasts, although some- 
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what degraded from the anechoic floor case. These results suggest that the phase- 
slope method described here may be less sensitive to the large random fluctuations 
that characterize single-point phase measurements under imperfect propagation condi- 
tions that are customarily encountered outdoors. 


CONCLUDING REMARKS 

A phase-slope method using two three-element microphone arrays has been used to 
predict the position of both pure— tone and broadband noise sources in an anechoic 
chamber. A methodology for quantitatively assessing both the accuracy and precision 
of such a source location technique has been introduced and illustrated. The effects 
of a change in ground impedance on the accuracy and precision with which a noise 
source can be located acoustically were investigated. When the phase-slope method was 
used, there was higher accuracy and more precision in estimates of source location 
with a small ground impedance than with a large ground impedance. Ground impedance 
effects were observed to be larger for low elevation angles than for high elevation 
angles. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
August 19, 1985 
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Array II 




Figure 2.- Microphone configuration. 



Array I 


Figure 3.- Array geometry viewed from top. 
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Frequency, Hz 

(a) Microphone pair B,A. 


Figure 9.- Phase spectrum corresponding to soft-floor m 

position 1 for broadband measurement 
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Figure 13.- Comparison of hard-floor and soft-floor measurements for array I 
obtained from pure-tone data and by theodolite. 
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(b) Elevation angle. Measured elevation angles for low source positions are 

0° for hard and soft floors. 


Figure 13.- Continued. 
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(a) Range. 

Figure 14.- Comparison of hard-floor and soft-floor measurements for array II 
obtained from pure-tone data and by theodolite. 










High, center 


(b) Elevation angle. 
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Measured elevation angles for low source positions are 
0° for hard floor. 


Figure 14.- Continued 
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Figure 15.- Comparison of hard-: 
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(b) Elevation angle. Measured elevation angles for low source positions are 
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Figure 15.- Continued 
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Figure 16.- Comparison of hard-floor and soft-floor measurements for array II 
obtained from broadband data and by theodolite. 
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(b) Elevation angle. Measured elevation angle for low, side source position 
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Figure 16.- Continued 
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(a) Pure-tone data. 

Figure 17*- Comparison of source positions measured from acoustic data with 

those obtained by theodolite. 
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Figure 17.- Concluded 
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(a) Source position 1. 


Figure 18.- Location accuracy and precision for four source positions. 






Figure 18.- Continued. 
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(c ) Source position 3. 


Figure 18.- Continued. 
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Figure 18 .- Concluded 
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